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ABSTRACT 

Using relativistic, steady, axisymmetric, ideal magnetohydrodynamics (MHD) we analyze the super-Alfvenic 
regime of a pulsar wind by means of solving the momentum equation along the flow as well as in the transfield 
direction. Employing a self-similar model, we demonstrate that ideal MHD can account for the full acceleration 
from high (^> 1) to low (<C 1) values of cr, the Poynting-to-matter energy flux ratio. The solutions also show a 
transition from a current-carrying to a return-current regime, partly satisfying the current-closure condition. We 
discuss the kind of the boundary conditions near the base of the ideal MHD regime that are necessary in order 
to have the required transition from high to low a in realistic distances, and argue that this is a likely case for an 
equatorial wind. Examining the MHD asymptotics in general, we extend the analysis of Hey vaerts & Norman 
and Chiueh, Li, & Begelman by including two new elements: classes of quasi-conical and parabolic field 
line shapes that do not preclude an efficient and much faster than logarithmic acceleration, and the transition 
a = a c after which the centrifugal forces (poloidal and azimuthal) are the dominant terms in the transfield 
force-balance equation. 

Subject headings: ISM: jets and outflows — MHD — methods: analytical — pulsars: individual (Crab Pulsar) 
— relativity — stars: winds, outflows 



1. INTRODUCTION 

The spin down luminosity of a rotating pulsar powers 
the emission of an associated synchrotron nebula by driving 
an initially Poynting flux-dominated wind (the Poynting-to- 
matter energy flux ratio is a 1) that flows along the open 
magnetic field lines. The density of the outflowing matter is 
enhanced by pair-creating cascades near the source and the re- 
sulting number density is sufficiently large to screen the elec- 
tric field along the flow and provide the necessary charge and 
current densities for ideal magnetohydrodynamics (MHD) to 
hold. 

At a distance r s s» 3 x 10 17 cm -specifically for the Crab - 
the ram pressure of the wind equals the total n ebula pressure 
and th e wind terminates in a standing shock ( Rees & Gunn 
fl97l . The "wisps" (which appear at distance ~ r s from 
the center) are an observational signature of the coupling be- 
tween the fast puls ar wind and the slowly expanding nebula 
llHesteretalJll995l) . By applying ideal MHD in the nebula, 
the matching shock conditions require that the pulsar wind has 
become completely matter-dominated at the position of the 
shock; the a function should be as low as a s = 0.003 just up- 
stream of the termination shock (Kennel & Coroniti 1984a,b), 
and the question arises as to how this transition from high to 
low values of a can happen, the so-called a problem. 

We note, however, th at the very exis tence of the shock 
has been questioned by Blandford (2002), who argues that 
no shock is necessary and that the observed emission is 
created by current dissipation in a Poynting flux-dominated 
wind. Furthermore, even if a shock indeed terminates the 
pulsar wind, abandoning the axisymmetry or the ideal MHD 
in the nebula results in undetermined cr, (Begermani ll998t 
Konigl & Granoll 120021) . Even so, a s could in general be 
small, and the interpretation of the wisps as ion driven com- 
pressions ( Arons 2002), a s well as t he fitting of the observed 
nebula emission (Bogovalov & Khangoulian 2002), provide 



additional support for the value a s = 0.003. 

The Crab-pulsar/nebula is the first and most analyzed object 
of this kind. Recent X-ray observations have shown that there 
exist other obje cts with common characteris tics [e.g., the Vela 
pulsar (Helfand. Gotth eTr&Halnernlb OOl ). the nebulae as- 
sociated wi th PSR B 1509-58 dGaensler et all200 j) . and PSR 
B 1957+20 ( Stapper s et al.l2003l) l. 

The study of approximately monopole magnetic field ge- 
ometries that are extremely inefficient accelerators gave the 
impression to the community that ideal MHD is in gen- 
eral unable to solve the cr problem, and many authors refer 
to the "cr parameter" and not to the "cr function". More- 
over, several_alternative (non-ideal MHD) attempts (see, 
e.g., lAronslfl998l for a review) gave unsatisfactory results, 
and the puzzle remains unsolved for almost three decades. 
However, there are in the literature r self-similar, ideal 
MHD solutions that start w ith a > 1 near the origin of 
the flow and reach doc « 1 (ILL Chiueh. & B egelman Jl992t 
IVlahakis & Konigl l|2003al). or even more efficient accelera- 
tors with oy, sa 0.5 (IVlahakis & Konigl l2003bT) . or « 0. 1 
( IVlahakis. Peng. & K5nTg]T2003l) . ("Actually, r self-similar are 
the only for the present known exact solutions of relativistic 
MHD in the sense that they satisfy the momentum equation 
in both directions: along the flow and in the transfield direc- 
tion.) We note in this connection that the main acceleration 
mechanism is not the centrifugal, but the magnetic. The for- 
mer may be important in the sub-Alfvenic regime of a disk- 
wind and results in poloidal velocity of the order of the initial 
Keplerian speed, while the latter - caused by the magnetic 
pressure-gradient force - results in much higher velocity that 
depends on the initial Poynting-to-mass flux ratio /j,c 2 . Em- 
ploying the cr function we may write = +cr oo), thus 
the maximum final Lorentz factor is attained in solutions with 
(Too = 0. In this paper we show that ideal MHD can account 
for the full acceleration from a ^> 1 to a « 0, the boundary 
condition 1/1 >> Afl at the base of the ideal MHD acceler- 
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ating regime is satisfied (here / is the poloidal current, A the 
poloidal magnetic flux function, and the field angular ve- 
locity). 

The transfield force-balance equation is the most impor- 
tant in determining the acceleration by means of controlling 
the field line shape and consequently how fast the quantity 
XS 2 B p decreases (here 03 is the cylindrical distance and B p 
the poloidal magnetic field). In the main part of the ac- 
celeration the transfield component of the electromagnetic 
force equals the poloidal centrifugal force that is propor- 
tional to the curvature of the poloidal field lines 1 /1Z. If 
this is still the case after the end of the acceleration - in the 
far-asymptotic regime where US/TZ = - the electromagnetic 
force e xactly vanishes yielding the solvability condition at 
infinity (Hevvaerts & Norman 1989; Chiueh, Li, & Begelman 
1199 11) . However, if during the acceleration phase sufficiently 
small a values are reached, the azimuthal centrifugal force 
- the other part of the inertial force - becomes important. In 
this regime it is the difference between poloidal and azimuthal 
centrifugal forces that equals the much smaller electromag- 
netic force in the transfield direction, and no solvability con- 
dition can be derived. 

We organize the paper as follows. In 5j2]we present the ideal 
MHD formalism focusing on the super-Alfvenic asymptotic 
regime where the puzzling transition from a ^> 1 to a -C 1 
happens. In *J3]we derive a novel z self-similar class of solu- 
tions describing the polar region of the wind and present the 
results of the integration. In |4] we analyze the asymptotic 
shape of the flow in relation to the acceleration, demonstrat- 
ing that the same results are expected in non-self-similar cases 
as well. In |5]we discuss the problem of the pulsar magne- 
tosphere, pointing out the problems that the force-free pic- 
ture faces, and suggesting qualitatively a scenario that seems 
plausible for creating an equatorial wind with the appropriate 
conditions for efficient acceleration from high to low-cr values 
satisfied. A summary follows in ^5] 

2. THE STEADY - AXISYMMETRIC MHD DESCRIPTION 

2.1. Governing equations 

The system of equations of special relativistic, steady, ideal 
MHD, consists of the Maxwell equations, Ohm's law, and the 
continuity, entropy and momentum equations: 

V-B = 0, VxE = 0, (la) 



V x B = 4nJ/c, V ■E = 4ttJ°/c, E = BxV/c, (lb) 
V-( A)7 ^) = 0, V-V(P/pl)=0, (lc) 

-7Po(^- V)(frV)-VP+(j°E+JxB) /c = 0. (Id) 

Here V is the velocity of the outflow, 7 the associated Lorentz 
factor, E , B the electric and magnetic fields as measured in 
the central object's frame, J°/c , J the charge and current den- 
sities, po >P the gas rest-mass density and pressure in the co- 
moving frame, while for the polytropic equation of state (with 
index V) the enthalpy-to-rest-mass ratio is 



£c 2 = c 2 + - 



r p 



(2) 



r-ipo 

Assuming axisymmetry [d/d(j> = 0, in spherical (r, #,</>) 
and cylindrical (z,t8,4>) coordinates with z along the rotation 



axis and the central object at (03 = ,z = z c )], five conserved 
quantities along the flow e xist. If A is the poloidal ma gnetic 
flux function, they are (e.g., Vlah akis & Konigll2003al) : 



the mass-to-magnetic flux ratio ^a(A) '■ 



47T7poVp 



the field angular velocity fi(A) = — - — — - , 

03 03 B p 

the specific angular momentum L( A) = £703 



the energy-to-mass flux ratio /i(A)c = £ 7 c 



O3f2B 



(3a) 
(3b) 
(3c) 
(3d) 
(3e) 



the adiabat Q(A) = — F , 
Po 

where subscripts p/<p denote poloidal/azimuthal components. 

An important combination of the above field line constants 
is the "Michel's magnetization parameter" 

AS1 2 

a M (A) = . (4) 

The physical quantities can be expressed as functions of A 
and the "Alfvenic" Mach number M = (^V p /B p )(4TTpoO l/2 ' 



Pa 



B = 



VA x 
03 

Yl = 

c 



AttM 2 ' 

= /i M 2 -(l-xj) 
7_ £ M 2 +x 2 -l ' 

x c 

M 2 B p V^ = ^7-M(l-4) 
c 



(5a) 
(5b) 
(5c) 
(5d) 



^7* A c c £ 7 x 

Here x = 03i7/c is the cylindrical distance in units of the light 
cylinder's lever arm on each field-streamline A =const, and 

1 /2 

xa = (LSl/pc 2 ) is its value at the Alfven point. 1 
Alternatively, using the equivalent to equation i5b\ . 

£7-^(1-4) 



M 2 -J^-x 2 — 

M"C7 P~il 

and the Poynting-to-matter energy flux ratio 



mQ,B d 



p-f7 



l + cr 



(6) 



(7) 



^Aiic 2 $7 

we may rewrite all quantities as functions of (A ,7) or (A, a). 

The two remaining equations are the Bernoulli and the 
transfield force-balance relations. The identity 7 2 - 1 = 

( 7 Vp/ c) + (jV<p /c) gives the Bernoulli equation 
7 2 -l = 

where 



£7~m(1-4) 
M"C7 " £(ii-&)x 2 

f7 2 03|VA| 



£7~M (1- 



(8) 



P^ac 3 



cr M 03 2 B ; , 



(9) 



1 We use the term "light cylinder" for the surface 03S1(A) = c, although in 
the general case where dQ/dA ^0 it is not a cylinder. 
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The projection of equation (licit along h = -VA/| VA| gives 
the transfield force-balance equation 



fcx + fix + fpx + /emi + fEMl + /em3 = , 



where: 



The azimuthal centrifugal term 



fcx = £,l Po-z-h-US 



Vl . fi 2 / My ^ x - 



C3 47TC3 V V p 



n-m. 



The rest of the inertial force density along h 



fix=-rPo&-[(V-V)V]-f C A 



.3 

Anil' 



-M 



is the poloidal centrifugal term. Here 1Z is the curvature 
radius of the poloidal field lines. 

The pressure gradient force density along h 
f P ±=-n-VP. 

The "electric field" force density 



1 



8ttC3 2 



«• V(U5 2 E 2 ) 



4ttTZ 



The "magnetic field" force density along ft 



1 



Tbx : 



8tt05 2 



T n- V(m 2 B 2 ) + 



4-ttTZ 4-irm 



-n-U5. 



The total electromagnetic force density in the transfield direc- 
tion if ei. + /b±) can be decomposed as 



1 



8ttG5 2 



«-V[05 2 (B 2 -E 2 )] + 



B 2 ,d-x 2 ) B 



ft, 



4ttK 

fEMl 



+ — '-n-® (10) 
4ttC5 



/eM3 



Altogether, they give the following form of the transfield 
force-balance equation: 



B 



4nTZ 



" (M 2 +x 2 -l) 



1 



47TG3 



8ttC3 2 
1 + 



n-V[05 2 (B 2 -E 2 )] + 

n-m-n-VP. (11) 



V, 



(Using equations (5), the latter becomes a second-order partial 
differential equation (PDE) for A.) 

At the origin of the flow (subscript i), xt <C Xa, Mf -C 1 — x\, 
and (using eq. (5b)) ji^il-xf) « p(\-x 2 A ). As the flow 
moves downstream, the functions x and M increase, and when 
the Alfven surface is reached they become x = Xa, M 2 = \-x\. 
The pulsar wind is still strongly magnetized in the neighbor- 
hood of the Alfven surface, and the field is close to force-free 
there: \f I± \ « \h p \\, or, M 2 < 1, implying x A < l. 2 The 
light cylinder x = 1 is located slightly after the Alfven surface. 
Next, the flow enters the super- Alfvenic regime. 

2 This is the case for trans-Alfvenic flows only. If the flow starts with 
super-Alfvenic velocity, xa could be in general larger than unity. 



2.2. The super-Alfvenic asymptotic regime 

In this regime x 2 ^> x\ w 1, M 2 ^> \-x\. Also thermal 
effects play no longer any role and the fluid is cold, £ w 1 . It is 
reasonable to assume that the above inequalities hold from the 
classical fast-magnetosound point downstream. At this point 
^ 2 - (B 2 -E 2 )/4ir p , and for cold flows, 7 ps [i 1 / 3 > 1 
;amenzindll986l) . 3 



(e.g.,Ei_ 

Employing the comoving magnetic field B 2 -E 2 
B\, the Bernoulli equation (|8j gives the exact result 



B 2 



£2= x 2 [x 2 -l + Af 2 (l-x 2 ) 2 /£ 2 ] B\ 
[x 2 -l+ M (l-x 2 )/£ 7 ] 2 ^ 



B 2 d-x z ) + 



(12) 



In the asymptotic regime where i> 1 the latter equation im- 
plies 



B 1 



1 2 



x 2 B 2 p _B\-{B 2 -E 2 ) l-l/ 7 2 



B 1 



B 2 (l-l/x 2 ) 1-1/x 2 



(13) 



Thus, the magnetic field is mainly azimuthal and the electric 
field approximately equals the magnetic field 4 



E=-|VA| =xB n 



-Ba 



1-1/7 2 



1/2 



1 - 1 /x 2 , 

The approximate form of the Bernoulli equation is 



(7 2 



1 



1/2 



: M7 
V-7 ' 



(14) 



(15) 



We are now in a position to examine the highly nonlinear 
transfield force-balance equation i ll 11 1. 

In that equation, the fs± and fg± are the dom inant 
terms that almost cancel ea ch other (e.g. IB o govalovl 1200 It 
Vlahakis & Konig fl l2003ah . The algebraic sum of their 
main contributions is the fEM\ term, which can be writ- 
ten as (using eq. fill ) -(l/87rC5 2 )n ■ V (C3B0/7) 2 . This 
term must be balanced by the other terms of equation 
\\ \\ . Ignoring fp±, we are left with the "poloidal curva- 
ture" term (B 2 /4t:TZ)(M 2 +x 2 - 1) = -fi±-fsM2 and the term 
-0B 2 /4^C5)[1 +(MV^/V p ) 2 ]n ■ 03 = -f c ±-f E m- Note that the 
last term, which henceforth we call "centrifugal", has con- 
tributions from the poloidal magnetic field (the part — /em3 
which is the dominant one in the c 3> 1 regime), and from the 
azimuthal centrifugal force —fc±. (The ratio of the two "cen- 
trifugal" parts is /ems/ fcx ~ M 2 /x 2 w cr, and the term fcx 
becomes important in the matter-dominated regime.) Thus, 
the asymptotic form of the transfield force-balance equation 



(M z +x z ) 



03 

K 



n ■ V 
203B 2 



®B 4 

7 



2„,2 



1 + 



M 2 7 



x 2 (7 2 -l) 



n-03. (16) 



This important equation, with the "cen trifugal" 
term omitted, was deriv e d bv | Chiueh etail ([991); 
Lvubarskv & E ichlerl (12001: lOkamotd J2002t) . while re- 
cently Tomimatsu & Takahashi (2003) derived the same 

3 Note that the flow is cold at the classical fast-magnetosound point when 
7:6 « M I/3 , or, (using eq. |5b| with Xj < x A , Mf < 1 -x\ ), 1 -x\ < 
(l-xf)fT 2 / 3 , implying x A ^ 1. 

4 In cases where E exactly equals -B^, eq . 1141 implies a "linear accel- 
erator" 7 = x 1 Contopoulos & Kazanas 2 002) . However, in the general case 
E does not exactly equals — Ba,, and eq. 1141 gives just the small difference 
between them. Note that t he E — —Ba, is not exactly satisfied in the numerical 
solution of Contopoulos, Kazanas. & Fendt 1 1999) either. 
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equation with the "centrifugal" term included, but the 
poloidal curvature term omitted. 

It is convenient to rewrite the Bernoulli and transfield force- 
balance equations in terms of the unknown functions (A , a) 
(also assuming 7 2 - 1 w j 1 ). The Bernoulli equation dl5> be- 
comes 5 7 = — £), or (using eqs. Q, and (9|). 



a er M 03|VA| 



ct m ra 2 ^ 



(17a) 



1 + cr /i A \ fi A 

The transfield force-b alanc e equation (I16> becomes (using 
eqs. E3, 00, and G7a|) 

C3 _ cr(l + er)G3VA ■ V Aficr c 2 03 • VA 

ft" i^VAj ^T"C3 2 2 |VA| ' ( } 

Equations dl7a> . and J17bl together with the equation for the 
curvature radius, 



K = 



IVAI 



V A - VA • V In |C3 VA I 



(17c) 



form a closed system for the functions a and A. A possi- 
ble solution of the system will describe the outflow down- 
stream from the classical fast-magnetosound surface (where 
the super- Alfvenic asymptotic conditions surely hold). 

Our task is to find a field geometry A(03,z), such that the 
flow accelerates smoothly from a ^ 1 (at the classical fast- 
magnetosound sur face) to a <C 1 (in the matter-dominated 
regime). Equation J17al shows that the key quantity 05 1 VA| = 
XS 2 B p should be a decreasing function of the distance from 
the source, along any field line. According to the model of 
Kennel & Coroniti (1984a,b) for the Crab nebula, the solution 
should give a s w 0.003, %^fi^ 10 6 , and B^ s « -2 x 10" 5 G, 
at the distance r s m 3 X 10 17 cm of the termination shock. At 
the classical fast-magnetosound surface 7/ « /i 1 / 3 = 10 2 and, 
from equation Q, <jj « 10 4 . 

2.2.1. General solution characteristics 

The current distribution 
The number density of the particles in the pulsar magneto- 
sphere highly exceeds the Goldreich-Julian value, and a tiny 
deviation from neutrality is enough to support the electromag- 
netic field. Assuming that in the upper (lower) hemisphere 
Z?0 < 0(> 0), the Poynting-flux is outflowing when the electric 
field points toward the rotation axis, corresponding to a nega- 
tive value of the charge density in the polar regions (J° < 0). 
For a highly relativistic poloidal motion, 



f> = — v • E = -— V(V x B) w — V„ ■ V X B w J» 

4-TT 4-TT 4-TT " 



(where a suffix "||" denotes the component of a vector along 
the poloidal field), and the flow that originates from the 
negatively-charged polar region is current-carrying (Ju < 0). 6 
The current-closure condition demands that each field line 
reaches a point where J» = 0, and after that enters the dis- 
tributed "return-current" regime /11 > 0. The poloidal current 

5 To a better approximation equation 1151 yields ^tS fi — 7 — [i/2-y 2 . 
The last term on the right-hand side is nonnegligible near the classical fast- 
magnetosound surface in the sense that the derivative of the latter equation 
(including the last term) along the flow gives (7~ 3 "A* -1 ) V ■ V7 = V ■ V£, 
implying 7/ = fi 1 / 3 , £/ = 1 -3/2/x 2 / 3 . However, fJ./2-y 2 <C A» and we neglect 
this term in our analysis. 

6 For the assumed quasi-neutral plasma, the sign of Jit depends on the tiny 
relative speed between the positively and negatively charged species. 



density is J p = ^ V x B^ = ^ V/ x 0, with 1 = JJ J p -dS = 
|n3Z?0, and the meridional current lines represent the loci of 
constant total poloidal current (/ = const). As a flow element 
crosses successive current lines, the absolute value of / de- 
creases and asymptotically vanishes {I^ = 0). 
An alternative scenario is that the current closes in a thin sheet 
at the end of the ideal MHD asymptotic regime, in which case 

/oo = /00(A) yo. 

The cr^l transition 
The field is force-free near the pulsar and the energy is Poynt- 
ing flux-dominated: a ^> 1. This remains the case near 
the classical fast-magnetosound surface, where a sa /i/7/ ~ 
/i 2 / 3 ^> 1 . As the Lorentz force accelerates the matter - mean- 
ing that the flow crosses current lines with decreasing |/| - the 
function a also decreases. The relation between a and / is (by 
combining equations I3dl . Q) 

0~ 2(—I) <7m 

= A-LJ±, (18) 

1 + cr Ail (1 

At a ~ a few the force-free approximation brakes down and 
the back-reaction of inertial forces on the magnetic field be- 
comes important. The final value of a could be in general 
very small (or exactly zero), depending on the small (or zero) 
value of loo ■ 

The cr > a c regime where /c±+/bm3 is negligible 
The "ce ntrifu gal" term in the transfield equation (the last term 
in eq. Ill7bl ) is negligible in the force-free regime. Since 
the flow velocity is mainly poloidal, it is expected that this 
term is also negligible inside the matter-dominated regime, in 
which case the poloidal curvature is controlled by the elec- 
tromagnetic force /em\ ■ Nevertheless, as the flow becomes 
more and more matter-dominated, the electromagnetic force 
/emi decreases (it vanishes at a = 0). As a result, a point 
a = a c is re ached where the two terms on the right-hand side 
of equation J17bl become comparable, and after that point the 
azimuthal centrifugal term takes over. 

For a > a c - when the centrifug al ter m is negligible - equa- 
tion dl7b> yields (using eqs. 0, 1 1 7 al ) 

'dcosiA _ ere 3 VA ■ V cr* A 

rfinoj ) A ~ |VA| 2 h~' ( } 

Here i? is the opening half-angle of the flow (the angle be- 
tween the poloidal magnetic field and the rotation axis). 
If the flow reaches a constant value (in principle different 
from field line to field line) a ^(A) such that > a c , then 
the right-hand side of equation dl9l is a function of A alone 
and c ost? would diverge logarithmically with 03 ( Oka motol 
120021) . The only way to avoid this divergence is to have 
c r oo(A)5 , A(A)/ri(A) =const. By using equations ( Bel , and 
the solvability condition at infinity 7oo(A)//oo(A) =const 
(Heyvaerts & Norman 1989; Chiueh et al. 1991) is recovered. 

The cr < a c regime where fc± is important 
At sufficiently small values of a (< cr c ), the curvature of the 
poloidal field-streamlines is controlled by the azimuthal cen- 
trifugal force 7 ; the electromagnetic for ce is much smaller than 
its inertial counterparts, and equation dl7b> gives 
03 cos $ 

TP"— (<0) - (20) 

In other words, in the a < a c regime we cannot neglect terms 
of order O (x^/x 2 ), O (03/7£.) in the transfield equation, sim- 
ply because the electromagnetic term is even smaller. 

7 The other part of the "centrifugal" force, namely the /em3> equals 
a fc± ^ fc± an d can be omitted in the matter-dominated regime. 
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In this regime the solution of the momentum equation is sim- 
ply V wconst and the motion is ballistic (in the poloidal plane 
the azimuthal and poloidal centrifugal forces cancel each 
other and the total inertial force vanishes). By decomposing 
a straight in three dimensions, constant velocity streamline in 
cylindrical coordinates, we get: V z = const, V < f > = c/x = c 2 /mft 
(from angular momentum conservation), and an increasing 

1 /2 

V m = (V 2 -V 2 -c 2 /x 2 ) . It can be shown that the poloidal 
streamline shape (the projection of the straight in three di- 
mensions streamline onto the poloidal plane) is hyperbolic, 



relations hold: 



V 2 -V 2 



mny (v 2 -v 2 ) 2 f (z-z )n 

c ) V 2 c 2 



with some constant zo- 

We emphasize that for flows with doo = <Too(A) < <j c , equa- 
tion d!9i should be replaced by equation J20i and there is no 
problem related to the divergence of cos In the absence of 
a solvability condition at infinity, any value CToo(A) < a c (and 
the corresponding current distribution, see eq. 1181 ) is a pos- 
sible solution. Even so, the most plausible case is Coo = cor- 
responding to a distributed return-current regime with 1^ = 
(Okamoto2002). 

The modified fast-magnetosound singular surface 
The classical fast-magnetosound surface is not singular when 
we solve simultaneously the Bernoulli and the transfield 
force-balance equations (it just separates the elliptic and hy- 
perbolic regimes of the PDE problem). The singular sur- 
face is the so-called modified fast-magnetosound surface (or 
fast-magnetosound separatrix surface) that lies inside the 
hyperbolic regime, coincides with a limiting characteristic, 
and plays the role of the event horizon for the pr opaga- 
tion of fast waves Q Tsinganos et al. 1996; Bogovalovl ll997t 
Vlahakis & Koniglll2003al) . The causality principle is satis- 
fied (disturbances in the asymptotic regime cannot influence 
the flow near the origin) only if the solution is trans-modified 
fast. Completely matter-dominated flows will have super- 
modified fast-magnetosound velocities (simply because the 
fast-magnetosound speed vanishes as a — > 0). 

3. THE Z SELF-SIMILAR MODEL 

Assuming that the poloidal field line shape is approximately 
cylindrical, i.e., B m <C B z \dA/dz\ <C dA/d(S, we may fur- 
ther simplify the system of equations dl7> to 



a cr M 05 dA 



1 + a fi A dm 1 

03 cr(l + cr)G3 d ( Aila\ c 
■ In ' ' 



K 



2 dm n V ctm J m 2 n 2 



(21a) 
(21b) 



To construct classes of analytical solutions, we shall make 
two assumptions: z is assumed a product of a function of A 
times a function of 03 (the z self-similar ansatz) 



Z = 03 $(A) C(BJ) 
(where 03o =const), and 

a = cr(03) . 



(22a) 



(22b) 



Following the algorithm described in 

IVlahakis & Tsinganosl (Q998), it is possible to separate 
the variables A and 03 in the system (121 1 only if the following 



A-A =--—Bom 2 <P F , 
(-t) 


(23a) 


/i = //()$. 


(23b) 


n=— $, 


(23c) 


fi (-F)A-Aq 
cm t A 


(23d) 



Here F,Bo,fio, and r are constants. The system J2 It be- 
comes 



= rG3- 

1 + cr (dU5 

m _m 2 (l + a)C F+2 md(c^c 
^z' 



2,2 



z 2 dm m L z 

Using the expression for the curvature radius 

<\_ ^v7 1+ b| x_3/2 



03 m 2 d 2 C { d( 



K z 2 ^ dm 2 \dm 



(24a) 
(24b) 

(25) 



we get finally a system of ordinary differential equations 
(ODEs) 

dm _ tC5(1 + (t) 

~d~c ~ 

l+cr 



da 



^T 2 m 2 1 + (T 



(26a) 



^C 3 

fi 2 T 2 m 2 



Fa 2 

c ' 



C03 2 



m 2 ( 2 <T 2 



(26b) 



The latter system can be easily integrated; the only difficulty 
is that the solution should cross a s ingu lar point that appears 
when the denominator in equation (I26bt vanishes. By rewrit- 
ing this denominator as 

[i\T 2 m 2 Hv^ 1 



m 2 ( 2 <7 2 



B 2 -E 2 



Anpoc 2 



(27) 



it is evident that the singular point corresponds to the mod- 
ified fast-magnetosound point, where the phase speed of the 
fast-magnetoso und waves propagating alo ng 03 is zero (see 
Appendix C in Vl ahaki s & Konigl 2003a). 

Using equation J22a> we may find VA and hence the mag- 
netic field, 



B = 



B m 2 ^ F ( a 



m * 

, , z+r— 03 . 
03 2 \l + a z J m l+cr 

whereas the current density component J» is given by 
cB A cr ( F + \ mda~ l 



B^m^ 



4tt l + cr 
The magnetic flux function is 

r(l + cr) 



dm 



A -An 



i-F). ^ 



(28) 



(29) 



(30) 



3.1. Results 
3.1.1. Scaling laws 
It is instructive to derive simple analytical scalings from 



equations (I26> : they turn out to be in very good agreement 
with the results of the numerical integration presented below. 
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Before analyzin g sep arately each regime, we rewrite schemat- 
ically equation J26bl as 

da _ -TZ l + B l -C 

dC K 2 -B 2 ■ { ' 

All the terms appearing on the right-hand side are positive (for 
a > t and F < 0). It is important to note that the terms 1Z\, 
IZ2 come from the poloidal curvature term of the transfield 
equation, the B\ and B2 have electromagnetic origin, while C 
comes from the "centrifugal" term. 

The sub-modified fast regime a ^> a m f 
In this regime the denominator in equation i31\ is negative, 
IZ2 < B2, i.e., the term Bi dominates. Since da/dC < 0, 
the numerator should be positive. Thus, the dominant term 
is the B\ (the only positive term). As we expected, in the sub- 
modified fast regime the magnetic field dominates over the 
inertial terms. One may write da/dC, « —Bi/Bj- The inte- 
gration gives a oc ( F and we finally have (by employing eq. 

Z oc $C5 1/T , a oc $ V . (32) 

Equation @3 gives 7 t | w Thus, for F < -1 the flow 

is current-carrying (Jn < 0) near the origin. 

The super-modified fast regime a <C <J m f 

In this regime, IZ2 > 62. Since da/dC, < 0, the inertial terms 
dominate in the numerator as well, resulting in da/dC, « 
-{IZi +C)/72-2- As long as a > a c we get that not only the 
poloidal curvature terms are important and should not be ne- 
glected, but also that the transfield force-balance equation re- 
duces to 03/7?. = 0. In other words, it is not correct to substi- 
tu te Q5/7£ sa (th e "pseudo-force-free" condition according 
to lOkamotol2002l) in the transfield equation and keep only the 
other terms. On the contrary, the XS/TZ w is the transfield 
equation itself. By including the "centrifugal" term we sim- 
ply modify the above conclusions and the transfield equation 
reduces to equation J20i . 

Thus, for a > a c the solution of the transfield equation is sim- 
ply a straight line C ~ (r s /®o) [l + (ct ! /t)(03-05 j )/05 j ] . By 
employing equation ( I26ai we finally get 

/ a, 05-05, \ 03 
z=r t $[l + — ), a = a s -. (33) 

V t 05., / m s z 

(In the most general, i.e., non-self-similar case, the solution is 
the conical la shape that we discuss in §0]) 
Equation <|29) gives J { \ w - c -^f^. Thus, for F > -2 the 
current density component J11 becomes positive at large dis- 
tances. The combination — 2 < F < — 1 corresponds to a tran- 
sition from current-carrying (Jh < 0, at a ^> 1) to a return- 
current regime (Jn > at a -C 1), required by the current- 
closure condition. 

The azimuthal centrifugal force becomes important after the 
point where the two terms on the right-hand side of equation 
d24bi become comparable. Using the scalings of equation 
( l33l . this happens at Tfi 2 x 2 a 2 , so a c « c^^fr /t5 m fQ,. 
If we continue the integration for sufficiently large z, the fast 
decrease of a stops at the point where a w r. At this point the 
product B p 6j 2 is of the order of (-F)A (see eq. 1301 ). As we 
discuss in §|4]the flow could continue to be accelerated after 
that point, but with a much lower (logarithmic) rate. 

3.1.2. Numerical integration 

Numerically we may integrate the equations J26t as fol- 
lows. Suppose that we examine the flow along the field- 
streamline A such that $(A)= 1, and \i = ^0 = 10 6 , Cl= 100 rad 



s" 1 (for which eq. 12 3d gives the value of 05o = c/57 = 3 x 10 8 
cm, and the independent variable of the ODEs is £ = zfi/c). 
We give the model parameters (r, F), and a trial value for the 
a function at the modified fast-magnetosound singular point, 
a = a m j. At thi s point, both the numerator and denominator of 
equation (I26bi vanish, and we are able to find the values Q m f 
and 05 m /. Using l'Hopital's rule we find the slope (da/dQ m j 
and start the integration from the singular point downstream. 
At some value of the independent variable £ = ( s we find 
<t = a s (= 0.003). At this point the distanc e sho uld be z = r s 
(= 3 x 10 17 cm), and according to equation ( 12 2 at the equality 
Q = r s /05o (= 10 9 ) should also hold. If the latter equality is not 
satisfied, we change the trial value a m f and repeat the whole 
procedure until we find Q s = 10 9 . For the correct er m y we are 
able to integrate downstream until a w a s (at the distance r s , 
just before the termination shock), as well as upstream until 
a = tjf w 10 4 , i.e., until the solution encounters the classical 
fast-magnetosound surface. (Using the value of the azimuthal 
magnetic field B$. s = -2 x 10~ 5 G and eq. 1281 . we find the 
constant Bo-) 

Solution a 

The representative solution a corresponds to the set of param- 
eters (t = 10~ 5 , F = -1 .9). For a m f = 0.9 the integration gives 
a s = 0.003 at r, = 3x 10 17 cm. Figure^shows the field lines as 
well as the current lines in the poloidal plane. The line shape 
is quasi-conical with small opening angle consistent with the 
model assumption B z ^> B m . The current density component 
/|| w /; changes sign from negative close to the classical fast- 
magnetosound surface to positive at larger distances, satisfy- 
ing (at least partl y) the current-closure condition (according 
to the analysis of SI3.1.1I this is expected for -2 < F < -1). 
Figure |2] shows the function a and the transition from 
Poynting- (a ~ 10 4 ) to matter- (a = 0.003) dominated regime. 
The a function scales as ~ z F = z~ 16 at a > a few (force- 
free regime) and as ~ 1 jz when the flow becomes matter- 
dominated. 

The components of the magnetic field are shown in figure [3] 
The main component is the azimuthal one, and at the distance 
r s it becomes 2 x 10~ 5 G. It is evident from the figure that the 
self-consistency condition B z ^> B m is everywhere satisfied. 
Figure|4]shows the transfield components of the various force 
densities as functions of z/r s along the reference field line 
$ = 1. As expected, the dominant forces are fg± ~ JnB^/c 
and /e± = J°E /c, and they almost cancel each other. For 
z/r s < 2 x 10~ 2 , fs± > and fg± < 0, corresponding to /|| < 
and J° < 0. For larger z/r s both fg± and fyx change sign and 
the flow enters the return-current regime J11 > 0, where the 
charge density is positive (J°/c > 0). 

The forces f IX +/em2, /emi , and f c ± +f E M3 ar e prop ortional 
to the first, second, and third terms in equation dl7b> . respec- 
tively. At small heights fj± + /emi < and the curvature ra- 
dius is 1Z > 0: its exact value is controlled by the electromag- 
netic forces \fj± + /emi] ~ \/emi \- However, at z/r s = 0.4, we 
encounter the point a = a c w 8 x 10~ 3 , where the "centrifugal" 
force —fc± -/em3 (that mainly consist of the azimuthal cen- 
trifugal part fc±) becomes important. After that point 7Z < 
and the transfield equation becomes \ fj± + /em2\ ~ —fcx- Th e 
small difference between these two terms is the electromag- 
netic term \f E m\- 
Solution b 

Figure [5] shows another solution similar to solution a. The 
only difference is that here r = 10~ 6 . Equations J26i show that 
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FIG. 1. — Solution a, corresponding to r = 10"" 5 and F = —1.9. (a) Three field lines (<J> = 0.25,1,4 - dashed lines) and four poloidal current lines 
[(-C5B^)i .2.3 .4 = 9 X 10" , 10 10 ,5 X 10' , 3 X 10'cgs - solid lines] for the solution a. As the flow moves along a particular field line it crosses current lines 
with de creas ing |/| = (c/2)(— KiB^), meaning that it is accelerated and the Poynting-to-matter energy flux ratio decreases. The field line shape is close to conical; 
see eq. 1331 . The current-carrying regime Ju < near the rotation axis, and the return-current regime /m > at larger cylindrical distances are shown, (b) The 
same field lines w ith a n enlarged scale, focusing on the region close to the classical fast-magnetosound surface. The field lines have 7Z > and the shape is 
parabolic; see eq. 1321 . A characteristic of the force-free regime is Bp \\ J p , e.g., the current line —UlB^ = (—GSB^i is initially tangent to the middle field line 
<J> = 1. (This is the case for z < 3 X 10 14 cm, for which fig.[2|shows that a > 9.) 



for different r, the value a m f remains roughly the same, while 
C3 m / scales as (JS m f)b = (®m/)a X r fl /r;,. Thus, the cylindrical 
distances in solution b are ten times larger compared to the 
ones for solution a. 
Solution c 

Figure [6] shows a solution corresponding to (r = 10~ 6 ,F = 
-2.5). In this case the constrain on a at r s implies o m j = 1.5. 
Because F < -2, the current density component Jh remains 
negative all the way from the source to infinity. If we con- 
tinue the integration to sufficiently large distances (not shown 
in fig. |6j the poloidal current becomes parallel to the flow and 
the acceleration stops; the function a then approaches a finite 
asymptotic value CToo(A). 

4. ASYMPTOTIC STRUCTURE OF EFFICIENT ACCELERATORS 

The outcome of the magnetic acceleration is either ctqc = 
Coo (A) or (Too = 0. Although the second possibility incorpo- 
rates the current-closure condition and seems more natural, 
the first one is also possible (it requires that the current closes 
in a thin sheet at the end of the ideal MHD asymptotic regime, 
but it cannot be ruled out, simply because all MHD winds 
practically terminate at some finite distance). 

10 4 
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10 1 
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FIG. 3. — The magnetic field components as functions of z, along the 
reference field line <E> = 1 of solution a. 



FIG. 2. — The Poynting-to-matter energy flux ratio a as a function of z = 
Cc/H, along the reference field line <E> = 1 of solution a. 
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FIG. 4. — Transfield components of the various force densities as functions 
of z/r s , along the reference field line $ = 1 of solution a. The transitions from 
J\i < to /ii > (at z/r s = 2 X 10~ 2 , where f B ± vanishes) as well as from 
1Z > to TZ < (at z/r s = 0.4, where fj± + /em2 vanishes) are shown. 

Besides the value of a oo, another important issue is the accel- 
eration rate. In this section we focus on the field line asymp- 
totic shape in the regime a < 1 and its relation to the accel- 
eration. (The analysis holds for nonrelativistic flows as well.) 
Note that the discussion of the present secti on - a s well as 
the characteristics of the flow described in § 12.2.11 - is gen- 
eral and not restricted to the self-similar solutions presented 
in § |3] It applies in the asymptotic regime a < 1 where the 
motion is mainly poloidal, <C V p . Equation d3bt yields 
-C3Z?0 oc U5 2 B p , or, (using eq. (HQ), a/(l + a) oc G5 2 B p (see 
also eq. I17al ). 8 Thus, an accelerating flow corresponds to a 
decreasing XB 2 B p , a function that depends on the line shape. In 
other words, how fast the acceleration takes place depends on 
the magnetic flux distribution, by the solution of the transfield 
equation. 

We already know that a field line shape of the form z = 
ci(A)C5''- (A) gives only logarithmic acceleration (Chiueh et al. 
1991). The deviatio ns from monopole magnetic field also re- 
sult in inefficient (Beskin, Kuznetsova, & Rafikov Hl998l) . or 
logarithmic (Lvubarskv & Eichler 2001) acceleration. (The 
monopole geometry itself - which however does not satisfy 
the transfield equation - gives ext remely inefficient acceler- 
ation: (Too = fi/joc = A* 2 / 3 > 1; lMiche]lH969l ) The loga- 
rithmic acceleration gives asymptotically er^ = 0, but over 
completely unrealistic (e xponentially large) distances. The 
numerical simulations by Bogovalov (2001) also give inef- 
ficient acceleration, possibly because the initial configuration 
is a split monopole field and does not change much during the 
simulation. 

The r self-similar soluti ons of Li et al. (1992) and 
IVlahakis & Konigll J2003al) give final equipartition = 
1, but the same model could give smaller as well 
JVlahakis & Konigl Il2003bt IVlahakis at al.ll2003l) . However, 
because it cannot capture the transition from positive to neg- 



8 In nonrelativistic flows, for V p approaching its maximum value V p , 
Vff,, a similar relation, namely -C5B^ rj U5 2 BpQ/V p0 o, holds. 



> 



ative curvature radius, it is impossible to obtain = 0. The 
asymptotics of the r self-similar model are straight poloidal 
field lines (05/7? = 0); as a result, the transfield component 
of the electromagnetic force asymptotically equals the "cen- 
trifugal" force (which could be very small, but definitely not 
exactly zero), and = a c . Even so, the value a c could be 
very small, and the model should be explored more carefully 
in that respect. 

In the z self-similar model presented in § [3] we kept the 
line shape free: z oc $(A)£(03), and the solution itself finds the 
shape through the function C(®)- The result is a quasi-conical 
shape, and the distribution of th e ma gnetic flux is different 
from the monopole field (see eq. I23al ). 

Collimation of relativistic outflows is in general much more 
difficult than in the nonrelativistic cases, firstly because the 
electric force almost cancel the transfield component of the 
magnetic force, and secondly because the effective matter in- 
ertia is larger (e.g., Bogovalov 2001). The result is negligi- 
ble curvature, 03/7?. « (e.g., Chiu eh! Li. & Begelmanl 1998). 
and only close to the origin where the Lor entz factor is sm aller 
than a few tens is collimation efficient ( Vlahaki s" & Konigll 
2003b). This result, together with the envisioned close to 
z = (53/ tan$(A) field line shape - we employ the term "type 
I conical" for this shape - led to the erroneous conclusion that 
acceleration is impossible when 03/7?. w 0. However, negli- 
gible curvature does not necessarily mean z oc 03; the most 
general case is a "type la conical" shape given by 



Z = Zo(A) + - 



03 



¥0. 



(34) 



tani?(A)' 

Differentiating equation J34> we get for the poloidal magnetic 
field 



, - cosi? £+sin?9 03 

WB P = 03 VA x 4> = — ; — . 

V /sini?-z sinu/G3 



(35) 



(Here primes denote derivative with respect to A.) It is evident 
from equation d35l that for z > 0, $ > the quantity 05 2 B p 
is a decreasing function of 03, meaning that the flow is ac- 
celerated. Magnetic flux conservation implies 2ir5A = B p 5S, 
where 6S is the cross-section area between two neighboring 
magnetic flux surfaces A, A + 6A (SS = 2ttG>5£±, see fig. Q. 
Since B p T8 2 / SA = 2irTB 2 /SS, it is clear that acceleration is pos- 
sible when - as the flow moves - SS increases faster than OS 2 . 
This is the case in type la conical flow, while in type I conical 
it is exactly SS oc OS 2 and the quantity B p 03 2 remains constant. 
The previous analysis shows that the argument of Chiueh et 
al. (1998) related to the global inefficiency of the magnetic 
acceleration in relativistic flows (based on the fact that the 
collimation is inefficient), can be circumvented. The conical 
la field lines are perfectly straight, satisfying the 03/7? = 0, 
nevertheless, they expand in the sense that SS increases faster 
than 03 2 . In other words, we may have expansion (increasing 
SS /03 2 ), without bending (03/7? j- 0). 

By combining equations fl35i and dl7aL we find that 
the (j function could change from a ~ 1 value at 03 ~ 
Az sin $/(A$ /sin?9 - 2a M /fS) (the conical la shape starts 
roughly at this distance), and becomes a = (T m ; n at distances 
03 > z sin 2 ?9/i?', with 



B p & 



o"m sini9 
~\x~a¥ 



(36) 



Since the factor sin^/Ai? is not <C 1 (this would require the 
opening of the field lines asymptotically to change rapidly 
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FIG. 5. — Solution b, corresponding to r = 10 6 , and F = —1.9. (a) Three field lines (# = 0.25 ,1,4- dashed) and three current lines (solid), (b)-(c) The 
Poynting-to-matter energy flux ratio a and the magnetic field components for the reference field line 4> = 1 . 




from line-to-line which is unlikely to happen) the only case 
that gives <7mi n < 1 is the \i <tm- This is an impor- 
tant requirement for the above mechanism to work, and it is 
equivalent to B p 6j 2 ^> A in the super- Alfvenic regime and as 
long as the a is 3> 1 [in this regime /i ss -03 flB ^/W^c 2 ~ 
BptSpQ/^AC 3 = (B p U5 2 /A)a M ]- Hence, as lChiueh et all i ll 9981) 
state, the poloidal field lines need to be bunched within a small 
solid angle before expansion takes place. Contrary to their 
claim, however, the expansion is not necessarily sudden near 
the light cylinder, but extended to jc ^> 1. In §[5] we discuss 
how such a field configuration may be created (see also fig. 

The general asymptotic expansion of the field line configu- 
ration - in both, relativistic and nonrelativistic cases - is 

z = zo(A) + c l (A)m C2{A) , (37) 

for which we get 



dU5 -c\/c\ -CjlnOJ-^/cOro" 6 ' 2 ' 



Chiueh et al. (1991) examined two types of parabolic cases: 
type I (zq = , C2 =/ 1 , c 2 = 0) when U5 2 B Z is a field line constant, 

and type II (z = , C2 =/ 1 , c 2 =/ 0) when the U5 2 B : slowly (log- 
arithmically) decreases and asymptotically vanishes. They 
also examined the conical type I case (z = 0,C2 = 1). Be- 
sides the conical type la case, new types of parabolic cases 
corresponding to z ^ should be added in the analysis. We 
henceforth call them type la, for C2 =const^ 1 , and type Ha, 
for C2 = q(A) ^ 1. These new parabolic types together with 
the type la conical, are much faster than logarithmic accel- 
erators in relativistic as well as in nonrelativistic flows. Ta- 
ble n summarizes the various possible geometries in the ac- 
celeration phase cr c < a < 1 (first column) combined with 
the possible final results of the acceleration - including pos- 
sible acceleration after the a = a c transition - and the far- 
asymptotic final geometries. The final current distribution 
and the dominant forces in the transfield direction are also 
shown. Compatible combinations of the geometry in the ac- 
celeration phase cr c < a < 1 and the far-asymptotic final state 
are marked with a "^/". For example, if the flow has a type 
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Table 1 

Possible asymptotic geometries and the corresponding acceleration^ 

The far-asymptotic regime (after the end of the main acceleration) 



a — > < (Too < a c 

1^0 I oo =I oo (A)^0 

flJ- = ~fcX iflX = -/cx)oo 

parabolic II or Ila hyperbolic 



The a c < cr < 1 acceleration phase 

film =-fi± 



c 00 =a c 
/ oo =/ oo (A)^0 

(Jem\ = -/cx)oc 
conical I or la 
or parabolic I or la 
or cylindrical 



Coo > OV 

3o =/oo(A) oc 7oo(A)* 

(jBMl)oo = 

conical I or la 
or parabolic I or la 
or cylindrical 



conical la: z = zo(A)+BJ/tan*i?(A) 
1 /CJ 2 B /; = ■&' / sin i? - zo sin tf/OS 



v 7 



v 7 



parabolic la: z = zo(A) + ci (A)C3 r2 , C2 =const 

1 /03 2 g; = (-cj /c i c 2 ) - (z'o [c i c 2 )0r C2 v/ 



v 7 



v 7 



v 7 



parabolic II: z = ci(A)05 C2<A) 
1/C!5 2 B, = (-c 1 /c l C2)-(c 2 /c 2 )lnU5 



V 



v 7 



parabolic Ila: z = zoG4) + ci(A)Of 2<A) 
l/WrB z = (-c , 1 /ciC2)-(c 2 /c2)lntn-(zo/cic 2 )Cr 



V 7 



* The table covers the asymptotics of nonrelativistic flows as well. 

* This is the solvability condition at infinity (Hevvaerts & Norman 1989; Chiueh et al. 1991). 



Ia conical shape (this is the most plausible case for relativis- 
tic MHD in which the field lines have negligible curvature 
GJ/72. <§; 1), the a function decreases according to equation 
dl7at and the shown in the Table^function XS 2 B p , and reaches 

a minimum value cr m i n at distances U5^> z sin 2 fl/d (see eq. 
1361 ). If the azimuthal centrifugal force is still negligible at 
these distances (i.e., er m j n > cr c ) the acceleration stops and the 
far-asymptotic regime is the one in the fifth column (most 
likely with unchanged conical Ia shape). If crmin ~ c c , then 
the azimuthal centrifugal force starts to play a role in deter- 
mining the poloidal curvature. One possibility is that the ac- 
celeration does not continue even with slightly different line 
shape, resulting in the fourth column case. The remaining 
case is that the a continues to decrease, and the line shape be- 
comes hyperbolic. Since the asymptotic form of a hyperbolic 
flow is again conical, the a function reaches a new minimum 
given from an equation similar to d36i . If the acceleration 
stops at this point, we get a final situation corresponding to the 
third column of Tabled However, even tiny deviations from 
an exactly radial flow Ui/TZ = could give a slow (logarith- 
mic) acceleration to the completely m atter-d ominated regime 
(Too = 0, as was argued by Okamoto (2002). This situation 
corresponds to a continuous acceleration under a parabolic II 
or Ila shape (second column). We note, however, that for re- 
alistic applications, the interaction of an outflow with its en- 
vironment brakes-down the ideal MHD conditions and this 
probably happens before the realization of this last part of the 
acceleration. 

Chiueh et al. (1991) examined the far-asymptotics of rel- 
ativistic MHD outflows, generalizing the w ork on nonrela- 
tivistic flows by Hevvaerts & Norman ( 1989). They assumed 
that the various quantities become functions of A alone (the 
characteristic of the far-asymptotic regime), neglected the 
poloidal centrifugal term 0{TB/1Z) as well as terms of order 
0{x\/x 2 ) in the transfield equation, and found that conical 
field lines satisfying (z/05)oo < oo should enclose a finite cur- 
rent loo = /qo(A ) oc 7oo(A) (the solvability conditi on at infin- 
ity, see 32.2. 1> . As already discussed in 32.2.11 this is the 



case when o-qo > a c , but for smaller values of Uoa the asymp- 
totic transfield equation reduces to the balance between the 
poloidal and azimuthal centrifugal forces [i.e., they are ex- 
actly the ©(x^/x 2 ), 0(63/7?.) terms that remain in the trans- 
field]; the electromagnetic term is much smaller, but never- 
theless not exactly zero and no solvability condition can be 
derived. Our self-similar solutions fall in this category: the 
condition Ioo{A) oc 7oo(A) is not satisfied, because < a c . 

As we see in Tabled the solvability condition at infinity 
is mandatory only for the fifth column. As a result we may 
have, e.g., nonrelativistic conical or type I/Ia parabolic flows 
that carry poloidal current 1^ = Ioo(A) (with negative or posi- 
tive dl oo jdA when they are in the current-carrying or return- 
current regime respectively), provided that < a c . 

As we already stated, the analysis of this section applies to 
nonrelativistic flows as well. A first step towards finding non- 
relativistic MHD solutions with Goo ~ is the supe r - mod - 
ified fast-magnetosound solution of Vlahakis et al. (2000), 
which ends slightly after the modified fast singular point with 
a = 0.005. The full acceleration (doo = 0) should be examined 
using a model that allows a transition from positive to nega- 
tive curvature radius (the asymptotic transfield force-balance 
equation for nonrelativistic flows implies that the curvature ra- 
dius has the sign of the —Jn, thus changing sign when the flow 
mov es from the curr ent-carrying to the return-current regime, 
e.g.. lOkamotol200l . 

5. THE PULSAR MAGNETOSPHERE/EQUATORIAL WIND 

The most efficient acceleration mechanism that we de- 
scribed is the one based on the conical Ia field line shape. The 
Z self-similar solutions are just a numerical example confirm- 
ing that the mechanism indeed works. Although the presented 
solutions describe an outflow very close to the polar axis, the 
same exactly physical mechanism is expected to work for the 
equatorial wind as well. Only practical difficulties do not al- 
low us to derive semianalytical solutions in that regime. 

As we discuss in § @J an important requirement for the 
above mechanism to work is that fj, 3> ctm- In fact, as we argue 
below, such a configuration is most likely to happen near the 
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equator. 

The pulsar magnetosphere has been traditionally consid- 
ered as being force-free, and inertial effects are usually ne- 
glected. On the other hand, it is generally accepted that it 
is because of the inertia that the magnetic field line topology 
becomes open and an azimuthal magnetic field component is 
created. Moreover, the pair creati on and their accelera tion 
to high Lorentz factors (e.g.. lDaughertv & Hardingll9 82) in- 
volve non-ideal MHD effects that may be important in the 
final form of the magnetosphere. A question arises: is the 
f orce-free a good approxim ation ? 

ICont opoulos et al. ( 199^) have solved the force-free prob- 
lem (a lthough wit h the c urrent-closure not fully satisfied) and 
lOgura & Koiimal J2003I) have repeated their numerical work 
with higher resolution, deriving similar results. The usual 
problem with the force-free dynamics is that the condition 
B 2 -E 2 > is not always satisfied, since the Bernoulli equa- 
tion that relates the electromagnetic field with the Lorentz fac- 
tor of the system of reference where the electric field van- 
ishes, is omitted (see eq. II 21 ). In regimes where B 2 < E 2 the 
drift velocity is larger than the light speed. This is indeed the 
case in the aforementio ned force-free solution. According to 
Ogura & Koiima (2003) the drift velocity becomes larger than 
the light speed at distance a few times the light cylinder (their 
figure 5 shows that typically this happens at cylindrical dis- 
tance « 3c/fi). Thus, the force-free assumption brakes-down 
much before the fast-magnetosound point (a = /U 2 ' 3 ). 

Another problem of the force-free solution is that it 
implies negative azimuthal velocities near the source. 
IContopoulos et all {1 999) demonstrated how we can find the 
flow speed in cases where the electromagnetic field is known. 
We need to solve the system of equations (coming from eq. 
(3b), and the second of eqs. 15 dl ) 



£7 1 



v ± = l _v ! ,j\i[_ 

cx c QBpti 2 ' 

-x—j = n(l-x 2 A ), 



(39a) 
(39b) 



together with the identity 1 - 1 / 7 2 = (V p /c) 2 + (V$ /c) 2 . The 
solution for 7 is actually equation Jl 2i . 9 The probl em re lated 
to negative azimuthal speeds comes from equation l !39at . The 
second part of the right-hand side is much smaller than unity 
very close to the source in a possible corotating regime, and 
should increase up to an asymptotic value equal to unity as we 
move to larger x (in order to have a negligible asymptotic az- 
imuthal speed). Practically speaking, the ratio 2\I\/VLB p US 2 
should be « 1 at distances x > a few. After that point 
(BpES 2 ) f ew ss 2|/|/fi, and this continues to be the case 
even in the superfast (non-force-free) regime where both the 
|/| and the B p XS 2 decrease. But, what happens relatively close 
to the light cylinder? The poloidal field lines follow the dipo- 
lar shape up to the point where they approach the equator, but 
then they remain parallel to it for some distance. We know 
that a relativistic flow is difficult to bend and this is a char- 
acteristic of a force-free regime as well. The curvature 1/7?. 
is of the order of (1/z) x O (max{l/x 2 , 1 /7 2 }), meaning that 
at distances x > a few, the poloidal field lines are practically 
straight. 10 As a result, the quantity XS 2 B p increases, because 

9 Since the back-reaction of the matter to the field is ignored, it is possible 
to find 1/7 = 0. This happens exactly in the positions were the drift velocity 
equals the light speed, i.e., where B 2 =E 2 . 

10 The z/TZ = (max{ l/x 2 ,! /■y 2 }) holds if we assume that h ■ V ~ l/z 



the surface between two neighboring field lines increases 
slower that cx 03 2 . For a force-free case where |/| remains con- 
stant, we get a decreasing ratio (V p /c)(2\I\/£lB p Gi 2 ). Since its 
asymptotic value should be unity, near the source is larger 
than unity, implying < 0. In fact Contopoulos et al. ( 1999) 
found negative azimuthal velocities. The only way to avoid 
this effect in the framework of force-free field, is to have very 
small initial poloidal velocity (V p <C c), meaning that our "ini- 
tial" surface is very close to the pulsar. However, in that case 
we should take into account non-ideal effects related to the 
pair creation region. 

A solution to this problem could be that the flow in the 
regime where the lines remain approximately parallel to the 
equator is decelerating, meaning that kinetic (or enthalpy) 
energy flux is transfered to Poynting, and the a function in- 
creases. In that case, the |/| = (c/2)G3|/?0 1 would increase and 
cancel the increase of B p US 2 in the expression of V^, resulting 
in V$ > 0. During the phase where the flow is approximately 
parallel to the equator, the B p XS 2 = 5A(2nXS 2 ) /6S increases al- 
most linearly with 03 (SS « 27r03fo with Sz ~ const). Since 
the 03|B^,| is forced to follow this increase (in order to be con- 
sistent with a negligible azimuthal motion), the matter energy 
part is decreasing linearly with 03, and an important (or even 
the most significant) part of the Poynting flux may be created 
there. As the flow decelerates, at some point it reaches 7^10. 
Only after that point the curvature of the poloidal field lines 
(1 /7Z) may be important and the lines start to expand to larger 
heights above the equator. The expansion in principle allows 
for acceleration (if SS increases faster than 03 2 ), and this time 
the electromagnetic energy is transfered to the matter. We 
consider this point as the origin of the ideal MHD regime that 
we examine in this paper. 

Another problem comes when we consider the "last" cur- 
rent loop (/ = 0) that closes on the equator. Along this equato- 
rial current line the azimuthal magnetic field vanishes. Since 
B 2 -E 2 = B\-B 2 p (x 2 -\), the inequality B 2 -E 2 > holds only 
in the unlikely case where the poloidal field also vanishes 
on the equator. The remaining case is to include non-ideal 
MHD effects (then the electric field is not equal to xB p and 
B 2 -E 2 = B 2 -E 2 could be positive). 

Although the details of the full magnetosphere remain to be 
explored, and this requires to solve the full problem (not only 
the non-force-free, but also the non-ideal MHD especially 
near the "last" current line that closes on the equator; not to 
mention the inner/outer gaps and the non-axisymmetry), the 
picture described above looks appealing and consistent with 
the observed Crab-pulsar equatorial wind. The scenario also 
explains why the quantities B p US 2 /A, 2\I\/A£l « h/<Jm, an< ^ 
1 /a are likely ^S> 1 just after the end of this "parallel to the 
equator" phase. These are the "initial" conditions for the ideal 
MHD phase. It may not be just a coincidence that the efficient 
and faster than logarithmic acceleration to Iow-ct values dur- 
ing the ideal MHD phase was possible for these "initial" con- 
ditions only. Note also that the part of the total open magnetic 
flux that contributes to the "parallel to the equator" phase, is 
~ 50% (for exactly dipolar field, 45% of the open magnetic 
flux correspond to field lines that have B z < when cross the 
light cylinder). The rest ~ 50% of the open field lines likely 
follows a dipolar shape till x ~ a few, and then becomes con- 
ical la (again 05/7?. is negligible). Since there is no obvious 

in eq. 1111 . Very close to the midplane, the field may change in a small scale. 
In addition, the pressure gradient term may be important. It is because of 
these two effects that the field lines do not cross the equator. 
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reason for an increasing B p U5 2 , the acceleration gives asymp- 
totically (Too ~ 1 or (ju if the initial a value is cr, 3> 1 or 07 -C 1, 
respectively. 1 1 

The "initial" conditions B p XS 2 A are in principle possi- 
ble in other astrophysical settings as well, such as gamma- 
ray burst sources or AGNs. For example, a corona or a ther- 
mal/radiation pressure field may keep the field lines that em- 
anate from a disk parallel to each other (at least for some dis- 
tance), in which case the product B p U} 2 again increases and 
a matter energ y-to-Poynting flux transformation is r ealized 
(e.g.. lVlahakis & Konigl l2003bllVlahakis at alJ200l . 

6. SUMMARY 

In the first part of the paper (§ l2.1ll2~2l we derived the ideal 
MHD equations that describe the super-Alfvenic asymptotic 
regime. The resulting system of equations remains intractable 
enough, and we could only find a z self-similar model (§ 
describing the regime near the rotation axis. We regard the 
solution a as the best representative result of that model. It 
includes all the expected characteristics of an efficient MHD 
accelerator, in particular: 

• The transition from a = 10 4 to a s» 0, showing that ideal 
MHD can account for the full acceleration, resulting in 
a matter-dominated flow. 

• The transition from a current-carrying to a return- 
current regime, showing how the poloidal current lines 
close. 

• The transition from positive to negative poloidal curva- 
ture radius. When the flo w beco mes sufficiently matter- 
dominated (cr < a c ; see 32.2. 1> it has a negative curva- 
ture radius and the line shape becomes hyperbolic. This 
is simply the result of angular momentum conservation: 
the azimuthal component of the velocity decreases, and 
the constancy of the total velocity implies an increas- 
ing cylindrical component. (This small deviation from 
a conical line shape could be enhanced in the termi- 
nation shock.) The most important implication of this 
effect is not the hyperbolic line shape itself (which still 
remains practically conical), but the new status quo in 
the transfield force-balance. It is now the poloidal and 
azimuthal centrifugal inertial forces that dominate and 
their difference is the much smaller transfield electro- 
magnetic force component. 

• The transition from sub- to super- modified fast- 
magnetosound flow. This surface should be crossed 
since in the matter-dominated regime the fast- 
magnetosound speed becomes progressively smaller (it 
vanishes at a = 0). In this sense it is a signature 
of an efficient acceleration. It is also related to the 
causality principle: only in the super- modified fast- 
magnetosound regime none of the MHD waves can 
propagate upstream and reach the origin of the flow. 

We tried to present our numerical solutions in accord with the 
model of Kennel & Coroniti for the Crab nebula. In solution 
a the transition from a = 10 4 to a s = 0.003 happens over less 

1 1 The creation of bunched poloidal field lines near the polar region cannot 
be completely ruled out. For example, the pressure of the created radiation 
may push the lines towards the axis. Another reason may be the Poynting 
flux (c/4tt)EuB^ associated with the parallel electric field inside the non- 
ideal MHD regime, which points toward the axis for 9 > 55° where Em < 0. 



than five decades in distance from the pulsar, and the value 
o s is re ached at distance of 3 x 10 17 cm. As we discussed in 
33.1.11 a decreases as oc l/z in the non-force-free regime a < 
a few (and even faster when a > a few). 

The assumption B~^> B m - required in order to derive the 
z self-similar model - does not allow us to examine the wind 
near the equator. Nevertheless, the presented numerical so- 
lutions serve as the first examples of low-er asymptotic flows, 
and, more importantly, they guided us to think which is the sit- 
uation in the most general (non-self-sim ilar) ca se. Our main 
conclusions (described and justified in § 12.2. H and are in 
summary: 

• Type la conical flows (eq. 1341 ) are efficient accelera- 
tors. The a decreases much faster than logarithmically. 

• The same for the parabolic la, Ha (they likely apply to 
nonrelativistic flows). 

• For sufficiently small asymptotic values of the a func- 
tion, a transition a = a c exist, where the azimuthal 
centrifugal force is comparable to the electromagnetic 
force in the transfield direction. The acceleration may 
or may not continue after this transition. In both cases, 
however, the solvability condition at infinity (Hey vaerts 
& Norman 1989; Chiueh et al. 1991) is not manda- 
tory. The latter should be satisfied only in cases where 
Coo > &c- Note that for the equatorial pulsar wind 
the ratio fEMi/fc± ~ ax 2 /j 2 is ^ 1 at the position 
of the shock (for a w 0.003, x ~ 10 9 , and 7 w 10 6 ). 
However, the solvability condition may not be satis- 
fied, if the value of the a required at the shock is not 
the asymptotic CToo, but the a function still decreases at 
the position of the shock. Moreover, in applications to 
jets related to other astrophysical settings, the transition 
er = er c may be important. 

• In flows where the acceleration continues to a < a c , the 
line shape is hyperbolic and the curvature radius of the 
poloidal field-streamlines becomes negative. 

A required condition for the above acceleration mechanism to 
work is [i om- As we discuss in §|4] this condition means 
that at the inner super-Alfvenic regime B p XS 2 should be ^> A, 
and the current that flows along the field lines should be much 
larger than the typical values obtained in monopole solutions, 
|/| 3> Ail. We argue in §|5]that this condition is likely satisfied 
for the equatorial wind. The scenario presented in that section 
is outlined below with the help of the figure The various 
transitions/regimes as we move along the reference (thicker) 
field line, marked with the numbered vertical dotted lines are: 

- Origin — > 1: This is the most uncertain regime, and very 
important problem for future research. We can only speculate 
what may happen there. The poloidal field possibly starts with 
a dipolar shape, however, inertial effects (associated with the 
azimuthal velocity) create an open magnetosphere. (Note that 
even under force-free conditions the last closed field line does 
not in principle intersect the equator at the light cylinder; see 
Uzdenskv 2003). Non-ideal MHD effects accelerate created 
pairs to 71 ^> 1. One has to solve the dynamics by includ- 
ing the back-reaction from the various emission mechanisms 
to yield the exact value of 71 . In principle, this value could 
correspond to matter-dominated flow ui < 1 (although the 
poloidal field is strong, the Poynting flux associated with the 
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EB<f, product may not be larger than the matter energy flux). 
The initial conditions after the pair acceleration in the non- 
ideal MHD regime may correspond even to super-Alfvenic 
conditions in which case the light cylinder has not a particular 
physical meaning. 12 

- 1 : Here x ~ a few. This may be the base of the ideal MHD 
regime (although there should be non-ideal effects near the 
midplane even after that point, see § [5jl. Downstream from 
this point the super-Alfvenic asymptotic analysis holds. Ac- 
cording to equation ( fT3t . xB p (— B^){\ - 1 /2"f 2 + l/lx 2 ) « 

and |/| = (c/2)®\B 4> \ w (fl/2)B p U5 2 . Also the poloidal 
curvature is G3/7£ -C 1 and the poloidal field lines are approx- 
imately straight and roughly parallel to the equator. 13 

~ 1 — ► 2: The product B p T8 2 increases as 2nTB 2 5A/5S = 
US6A/Sl±, where Sl± is the distance between two neighboring 
field lines A, A + 8A (see fig. 0. Since the (SE oc B p U5 2 in- 
creases, |/| = (c/2)W\B^,\ should also increase in order to keep 
the inequality B 2 > E 2 true. Thus, matter kinetic energy (and 
possibly thermal energy) is transfered to Poynting. Figure 
shows a current loop that crosses the reference field line be- 
tween points 1 and 2. Here J± < and the associated Lorentz 
force (c/4ir)J±B c f ) decelerates the matter, increasing the a 
function. Also J\\ < 0, J° < 0, but the transfield component 
of the Lorentz force [fa±+fE± ~ (c/4n)(J°E + J\\B^)] re- 
mains small, result ing in Ui/TZ = O (max{l/x 2 , 1/7 2 }), (e.g. 
iChiuehet alii 9981 

- 2: Here the Lorentz factor is sufficiently small (< 10) for 
the field lines to start slightly bending. This is the base for 
the accelerating ideal MHD regime described in the previ- 
ous sections. This point corresponds to an "initially" Poynt- 
ing flux-dominated wind. The "initial" conditions satisfy 

12 In cases where the flow starts with super-Alfvenic speed, the f2 is 
no longer related to the m atter rotati on. It is just the ratio cE/U5B p (e.g., 
Conto DouloJ1995tlVlahakis & Konigll2003bl) . 

13 Two other factors that could help the squeezing of the field lines in a 
regime close to the equator are gravity and ram pressure of infalling mate- 
rial from the environment of the pulsar. However, both effects are probably 
insignificant. 



B p G5 2 = (B p G5 2 )i » A [1 » a M and the a function is >• 1 . 
-2^4: The field lines have a parabolic shape: 5l±/{fi5A 
slightly increases, resulting in a tiny decrease i n B P US 2 . The 
a function decreases according to the equation (I17a> . For the 
poloidal current density, Jn < 0, J± > 0. 

- 3: Classical fast-magnetosound point, 73 m /i 1 / 3 , 0-3 w /i 2 / 3 . 

- 4: The a function is unity. Equation d!7al implies that 

B p m 2 = {B p m 2 )i/2. 

- 4 — > 6: This is the conical la phase (the discussion of §|4] 
refers to this regime). The field lines become z = Zo(A) + 
05/ tant?(A) (tangent to the dot-dashed lines in fig. 0. It is 
seen in figure how the conditions z > 0, •& > are real- 
ized. The B p XS 2 decreases according to equation (1351 . As a 
result, the a function decreases from w 1 to cr m i n given by 
equation d36i . 

- 5: The current loop that the reference field line crosses at 
this point has Jh = 0, i.e., here the current-carrying regime 
(7|| < 0) ends and the return-current regime begins (Ju > 0). 

- 6: Here the B p XS 2 /A is of the order of unity. The fast ac- 
celeration cannot continue after that point. However, due to a 
tiny deviation from a XS 1 11 = it may continue logarithmically 
(see lChiueh et al.ll9 98 and references therein). 

If at point 6 the value of erg is <j s = 0.003 while at point 
2 the flow is Poynting dominated (02 ^1), then the ratio 
(B p G5 2 )(,/(B p G5 2 )2 = 06 and the required value for the B p G5 2 /A 
at point 2 is l/<76 ~ 300. Assuming further that ideal MHD 
holds during the 1 — > 2 phase and that the flow started at 1 with 
a x < 1, we get (B p U} 2 ) 2 /(B p lS 2 )i = (fl*tB) 2 /(fl*tB)i = 1/tn. 
Thus, {B p TB 2 )\/A = ai/a^ w 300cti. For a dipolar field before 
the point 1, (B p XS 2 )\/A = 2, and we get a required value for 
a 1 = 2(76 ~ 0.006. However, it is not clear that all the above 
assumptions hold, so we cannot be sure for the derived value 
01. 

In general, the fast (faster than logarithmic) acceleration in 
magnetized outflows gives asymptotically <r m i n ~ A/B p XS 2 = 
cm/V' where B p XS 2 is the "initial" value at the Poynting dom- 
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inated regime (point 2 in fig. [7}. Besides the application to 
pulsar winds on which we focussed in this paper, the formula- 
tion presented is quite general and AGN or GRB relativistic as 
well as YSO nonrelativistic jet asymptotics can also be con- 
sidered (the asymptotic analysis presented in §|2]applies in all 
these cases, no matter if asymptotically a <C 1 or ~ 1, or even 
> 1). In addition, the z self-similar model of §|3] could be 
used for examining AGN or GRB jets in their super- Alfvenic 
regime. 
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